Gravitational Condensate Stars: An Alternative to Black Holes 
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A new solution for the endpoint of gravitational collapse is proposed. By extending the concept of 
Bose-Einstein condensation to gravitational systems, a cold, compact object with an interior de Sitter 
condensate phase and an exterior Schwarzschild geometry of arbitrary total mass M is constructed. 
These are separated by a phase boundary with a small but finite thickness I of fluid with eq. of 
state p — +pc 2 , replacing both the Schwarzschild and de Sitter classical horizons. The new solution 
has no singularities, no event horizons, and a global time. Its entropy is maximized under small 
fluctuations and is given by the standard hydrodynamic entropy of the thin shell, which is of order 
k B £Mc/h, instead of the Bekenstein-Hawking entropy, S BH = 4irk B GM 2 /he. Unlike black holes, a 
collapsed star of this kind is thermodynamically stable and has no information paradox. 
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Introduction. The vacuum Einstein eqs. of classical 
general relativity possess a well-known solution for an 
isolated mass M, with the static, spherically symmetric 
line element, 



ds 2 = -f{r)dt 2 + 



dr 2 
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{de 2 



sm 
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where the functions f(r) and h(r) are given by 



f(r) = h(r) = 1 



2GM 



(1) 



(2) 



in units where c = 1. The dynamical singularity of this 
Schwarzschild metric at r = with its infinite tidal forces 
clearly signals a breakdown of the vacuum Einstein eqs. 
The kinematical singularity at the Schwarzschild radius 
R s = 2GM is of a different sort, corresponding to an 
infinite blue shift of the frequency of an infalling light 
wave with respect to its frequency far from the black 
hole. Since the curvature tensor is finite at r = i? s , the 
singularity of the metric there can be removed by 

a suitable change of coordinates in the classical theory 
A classical point test particle freely falling through the 
event horizon experiences nothing special at r = R s . 

The physics at the event horizon may be quite different 
when h 7^ 0. Consider the same photon of finite asymp- 
totic frequency u> far from the black hole. Even if Tiuj is 
arbitrarily small, the infrared photon acquires a local en- 
ergy E = htuf~2, which diverges at r = R s . Hence there 
is no a priori small parameter controlling deformation of 
the local geometry near r = R s due to quantum effects. 

In the semi-classical approximation, when a massless 
field such as that of the photon is quantized in the fixed 
Schwarzschild background, one finds that the black hole 
radiates these quanta with a thermal spectrum at the 
asymptotic Hawking temperature, T H = h/SnksGM |Q. 
It is usually assumed that the backreaction effect of this 
radiation on the classical geometry must be quite small. 



However, detailed calculations of the energy-momentum 
of the radiation show that its (T\) and (T r r ) components 
have an / _1 infinite blue shift factor at the horizon [p|. 
The wavelengths contributing to these quantum stresses 
are of order R s and hence are non-local on the scale of 
the hole. One has only to insert such a diverging en- 
ergy density and pressure into the semi-classical Einstein 
eqs. (||) below to see that the geometry near r = R s 
is changed significantly from the classical Schwarzschild 
form. Unlike the classical kinematic singularity in (Q), 
such non-local semi-classical backreaction effects near R s 
cannot be removed by a local coordinate transformation. 

Furthermore, the inverse dependence of T H on M im- 
plies that a black hole in thermal equilibrium with its 
own Hawking radiation has negative specific heat and 
therefore is unstable to thermodynamic fluctuations [0. 
Energy conservation plus a thermal radiation spectrum 
also imply that the black hole has an enormous entropy, 
S BH ~ 10 77 k B (M/M Q ) 2 (|, far in excess of a typical stel- 
lar progenitor. The associated information paradox has 
been conjectured to require an alteration in the principles 
of relativity, or quantum mechanics, or both. 

In lieu of such fundamental revision of the laws of 
physics, it is reasonable to examine alternatives to the 
strictly classical view of the event horizon as a harm- 
less kinematic singularity, when /i / and the quan- 
tum wavelike properties of matter are taken into account. 
In earlier investigations which attempted to include the 
backreaction of the Hawking radiation in a self-consistent 
way, the entropy arises entirely from the radiation fluid 
§|J. In fact, 5 = 4 j^f T S BH , for a fluid with the eq. of 
state, p = np, becoming equal to the Bekenstein-Hawking 
entropy for k = 1. Despite such intriguing features, these 
models have huge (Planckian) energy densities near R s 
and a negative mass singularity at r = 0, so that the 
Einstein eqs. are not reliable in the interior region. 

On the other hand, instead of abandoning quantum 
mechanics, demanding that its principles apply to self- 
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gravitating systems implies that their specific heat oc 
((AE) 2 ) must be positive. This leads to the proposal that 
quantum black hole entropy is carried not by Planckian 
but by soft quanta with E ~ h/GM 0. More recently 
another proposal for incorporating quantum effects has 
been made, viz. that the horizon becomes a critical sur- 
face of a phase transition in the quantum theory, sup- 
ported by an interior region with eq. of state, p = —p < 
Such a vacuum eq. of state, first proposed by Gliner 
m for the cndpoint of gravitational collapse, is equivalent 
to a positive cosmological constant in Einstein's eqs., and 
does not satisfy the energy condition p + 3p > needed 
to prove the classical singularity theorems. 

In this Letter we show that an explicit static solution 
of Einstein's eqs. taking quantum considerations into 
account exists, with the critical surface of ref. Q replaced 
by a thin shell of ultra-relativistic fluid of soft quanta 
obeying p — p. Such a solution, lacking a singularity 
and an horizon is significant because it provides a stable 
alternative to black holes as the endpoint of gravitational 
collapse, possibly with different observational signatures. 

The assumption required for this solution to exist is 
that gravity undergoes a vacuum rearrangement phase 
transition in the vicinity of r = R s . Specifically, in this 
region quantum fluctuations on the scale R s dominate 
the T* t ~ T r r components of the stress tensor, which 
grow so large that the eq. of state approaches the most 
extreme one allowed by causality, p — p. As this causal 
limit is reached, the interior spacetime becomes unstable 
to the formation of a new kind of gravitational Bosc- 
Einstein condensate (GBEC) described by a non-zero 
macroscopic order parameter \ff. If |W| 2 is a constant 
scalar, it must couple to Einstein's eqs. in the same way 
as a cosmological term, and the eq. of state of the inte- 
rior region must be p v = — p v — V(\^\ 2 ). A suggestion 
for the effective theory incorporating the effects of quan- 
tum anomalies that could give rise to both this interior 
GBEC phase and the p = p shell has been presented else- 
where 10 . Here wc forego any discussion of the details 
of the phase transition or collapse process and consider 
only the compact, stable endpoint of gravitational col- 
lapse, by solving the static Einstein's eqs. with the spec- 
ified eqs. of state. For a recent review of investigations 
of other non-singular quasi-black-hole (QBH) models see 
ref. 0. 

Solution of Eqs. The eqs. to be solved are the Einstein 
eqs. for a perfect fluid at rest in the coordinates ([!]), viz. 



which ensures that the other components of the Einstein 
eqs. are satisfied. These three first order eqs. for the four 
unknown functions, /, h, p, and p become closed when an 
eq. of state for the fluid, relating p and p is specified. 
Because of the considerations above we allow for three 
different regions with the three different eqs. of state, 



I. Interior : < r < r\ , p = — p , 

II. Shell : ri < r < T2 , p = +p , 

III. Exterior : r% < r , p = p = . 



(5) 



At the interfaces r = r% and r = r 2 , we require the 
metric coefficients r, / and h to be continuous, although 
the first derivatives of /, h and p must be discontinuous 
from the first order eqs. (^) and (|4|). 

In the interior region p = —p is a constant from (Q). 
Let us call this constant p v = 3Hq/8itG. If we require 
that the origin is free of any mass singularity then the in- 
terior is determined to be a region of de Sitter spacetime 
in static coordinates, i.e. 



I. 



fir) = Ch(r) = C(1-H 2 r 2 ), < r < r x . (6) 



The unique solution in the exterior vacuum region which 
approaches flat spacetime as r — > oo is a region of 
Schwarzschild spacetime (0), viz. 



III. 



f(r) = h(r) = 1 - 



2GM 



Ti < r . 



(7) 



The integration constants C, Hq and M are arbitrary. 

The only non- vacuum region is region II. Let us define 
the dimensionless variable w by w = 8nGr 2 p, so that eqs. 
(§)-(|) with p = p may be recast in the form, 



dr 
r 

dh 
~h~ 



dh 



1 — w — h 

1 — w — h dw 
1 + w — 3h w 



(8a) 
(8b) 



together with pf oc wf jr 2 a constant. The first of eqs. 
(||) is equivalent to the definition of the (rescaled) Tolman 
mass function by h = 1 — p,/r and dp{r) = 2 G dm(r) = 
8itG pr 2 dr — w dr within the shell. The second eq. ( |Sb| ) 
can be solved only numerically in general. However, it is 
possible to obtain an analytic solution in the thin shell 
limit, < h <C I, for in this limit we can set h to zero 



on the right side of (8b) to leading order, and integrate 



it immediately to obtain 



-G\ = 4r -r [r (1 - h)] = -8nG T\ = 8ttG p , (3a) 
r z dr 

G\ = A.Z + L (h - 1) = 8nGT r r = 8ttG P , (3b) 
r) dr r z 



together with the conservation eq., 

Va r dr + 2/ dr ~ U ' 



(4) 



p (l + w) 2 

— ~ e < 1 

r w 



(9) 



in region II, where e is an integration constant. Because 
of the condition h <C 1 we require e <C 1, with w of order 
unity. Making use of eqs. (||) and (||) we have 



dr 
r 



-e dw 



(1 



(10) 
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Because of the approximation e •< 1, the radius r 
hardly changes within region II, and dr is of order 
edw. The final unknown function / is given by / = 
(r I fx) 2 (w\ I 'w) f (n) ~ (wi/w)f(ri) for small e. 

Now requiring continuity of the metric coefficients / 
and h at both r\ and r 2 gives the conditions, 



h( n ) = 1-H£ri ~e- 

2GM (l + w 2 ) 2 



% 2 ) 
/(ra) 



1 



1 
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Wl /(^l) 

u; 2 h(r 2 ) 



w-2 



C 



l + Wx 

1 + w 2 



(11a) 
(lib) 
(11c) 



Together with r%/ri from the solution of (|To|) this gives 
four relations among the eight integration constants 
(ri, r 2 , toi, W2, Hq, M, G, e). Hence the first four can be 
eliminated in favor of 7?o, M, C and e <C 1, and we have 
a four parameter family of static solutions. Assuming 
that (7*1, r 2 , Wi, W2) remain finite as e — ► 0, i.e. are of 
order e°, we obtain from (|o|) and ( |i"T| ) that 7*1 — -ff^ 1 — 
r 2 ~ 2GM , but 1 - ff n , r 2 - 2GM, wi - w 2 and C - 1 
are all of order e, while r 2 — ri is 0(e 2 ). 

Principal Features. If e > both / and h are of order e 
and approximately constant in region II, but are nowhere 
vanishing. Hence there is no event horizon, and t is a 
global time. A photon experiences a very large, 0(e~i) 
but finite blue shift in falling into the shell from infinity. 
The proper thickness of the shell, 



t = 



dr h 



dw \ 



(12) 



is small compared to R s . Although £ is arbitrary here, 
presumably it is fixed by microscopic physics and is not 
greater than the Planck length by more than a few or- 
ders of magnitude. If we assume that £ is a fixed multi- 
ple of the Planck length and independent of M, then 
e ~ {£/GM)i ~ {Mpi/M)i ~ 1CT 25±1 for a solar 
mass object, which certainly justifies the small e approx- 
imation in this case. With £ fixed, H^ 1 ~ 2GM and 
(7=1 + 0(e), the only remaining free parameter is M. 

The entropy of the shell is obtained from the eq. of 
state, p = p = (a 2 /8irG)(k B T/fr) 2 , where we have in- 
troduced G for dimensional reasons so that a 2 is a di- 
mensionless constant. By the standard thermodynamic 
relation, Ts = p + p for a relativistic fluid with zero 
chemical potential, and the local specific entropy density 
s{r) = a 2 k 2 B T(r)/4Trh 2 G = a(k B /h)(p/2irG)i for local 
temperature T(r). Thus s = (ak B /AirhGr) and the 
entropy of the fluid within the shell is 



S = 47r / sr 2 drh 2 



ak B R 2 1 . (w x 
nG \W2 



(13) 



Using W1/1V2 = 1 + 0(e) and ( |l2] ) this is of order 



S 



ak B Y R s e '' 



ak r 



M£ 



«5 K 



(14) 



Since the interior region I has p v — —p v , Ts vanishes 
there, as could be anticipated for a GBEC described by 
a macroscopic single quantum state. The entropy of the 
entire quasi-black hole (QBH) is given then by ( filf ) or 
(|l4|), which is of order 10 38 fc B (l/Lpi) for a solar mass 
object. This is some 38 orders of magnitude smaller than 
the Bekenstein-Hawking entropy for the same mass M, 
and 20±1 orders of magnitude lower than a typical stellar 
progenitor, which have entropies in the range of 10 57 fc B 
to 10 59 fc B for M/m N - 10 57 nucleons. Since w is of order 
unity in the shell while r ~ R s , the local temperature of 
the fluid within the shell is of order T H ~ Ti/k B GM, 
i.e. the typical quanta are soft. Because of the absence 
of an event horizon, the QBH does not emit Hawking 
radiation. In fact, the thermal wavelength of the soft 
quanta is of order the linear dimension R s of the object, 
so that quantum zero point fluctuations and finite size 
effects are competitive with classical thermodynamics. If 
w is somewhat smaller than unity, the shell need not 
emit any thermal radiation at all. A discussion of the 
corrections to the classical hydrodynamic relations this 
entails must await a fuller quantum treatment. With no 
thermal emission the GBEC remnant then would be both 
ultracold and completely dark. 

The extremely cold radiation fluid in the shell is con- 
fined to region II by the surface tensions at the time- 
like interfaces n and r 2 . These arise from the pres- 
sure discontinuities, Api ~ H 2 (3 + wi)/87rG and Ap 2 ~ 
— w 2 /327rG 3 Af 2 , and are calculable by the Israel junction 
conditions (ljjl^]. We find the non-zero angular compo- 
nents of the surface tension to leading order in e, 
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32ttG 2 M (1 + toi) 
1 w 2 



32irG 2 M (1 



(3 + iUi) ^wiji 

(t)' 



W2I 



(15a) 
(15b) 



at n and r 2 respectively. The signs correspond to the in- 
ner surface at r\ exerting an outward force and the outer 
surface at r 2 exerting an inward force, i.e. both surface 
tensions exert a confining pressure on the shell region II. 
The time component a l t — 0, corresponding to vanish- 
ing contribution to the Tolman mass function m(r) at 
the two interfaces. Since e _2 ^ [M /M p i)s , the surface 
tensions ([if]) are of order M~i and far from Planckian. 
Hence the matching of the metric at the phase interfaces 
r\ and r 2 , analogous to that across stationary shocks in 
hydrodynamics, should be reliable. Resolving the inter- 
faces will require going beyond Einstein's eqs. to a more 
microscopic description of the quantum phase transition. 
The energy within the shell (as measured at infinity), 



Ey 



4-7T 



pr dr ~ eM 



dw 



(l + w) ~ e 2 Af (16) 
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is also extremely small. Hence essentially all the mass 
of the object comes from the energy density of the vac- 
uum condensate in the interior, even though the shell is 
responsible for all of its entropy. 

Stability. In order to be a physically realizable end- 
point of gravitational collapse, any quasi-black hole can- 
didate must be stable H. Since only the region II is 
non-vacuum, with a 'normal' fluid and a positive heat 
capacity, it is clear that the solution is thermodynami- 
cally stable. The most direct way to demonstrate this 
stability is to work in the microcanonical ensemble with 
fixed total M, and show that the entropy functional, 



S : 



TiG 



rdr 



dr ) 



(17) 



is maximized under all variations of p(r) in region II with 
the endpoints (ri,^), or equivalently (id 1 ,w 2 ) fixed. 

The first variation of this functional with the endpoints 
T\ and r-i fixed vanishes, i.e. SS — by the Einstein eqs. 
(||) for a static, spherically symmetric star. Thus any 
solution of eqs. is guaranteed to be an extremum 

of S H] . This is also consistent with regarding Einstein's 
eqs. as a form of hydrodynamics, strictly valid only for 
the long wavelength, gapless excitations in gravity. 

The second variation of (^) is 



S 2 S 



ak B 


- / r 
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dr 



( 

\ dr 



(5^f dn ( + d^ 
r 2 h 2 dr V dr , 



(18) 



when evaluated on the solution. Associated with this 
quadratic form in Sp is a second order linear differential 
operator C of the Sturm-Liouville type, viz. 




dr 



dp\ 2 L dp 
dr J \ dr 



(19) 



This operator possesses two solutions satisfying C\ a = 0, 
obtained by variation of the classical solution, /x(r; r% , r2) 
with respect to the parameters {r\,T2)- Since these cor- 
respond to varying the positions of the interfaces, % does 
not vanish at (ri,r2) and neither function is a true zero 
mode. For example, it is easily verified that one solution 
is x a — ^ — w. However, we may set Sp = x "0; where 
ij) does vanish at the endpoints and insert this into the 
second variation (Esj). Integrating by parts, using the 
vanishing of Sp at the endpoints and Cx = gives 



S 2 S=- 



ak B 
4hG 



r dr [ — — 

dr 



dtp 
dr 



< 0. 



(20) 



Thus the entropy of the solution is maximized with re- 
spect to radial variations that vanish at the endpoints, 
i.e. those with fixed total energy. Since deformations 
with non-zero angular momentum decrease the entropy 
even further, stability under radial variations is sufficient 
to demonstrate that the solution is stable to all small per- 
turbations. In the context of a hydrodynamic treatment, 
thermodynamic stability is also a necessary and sufficient 
condition for the dynamical stability of a static, spheri- 
cally symmetric solution of Einstein's eqs. H]. 

Conclusions. A compact, non-singular solution of Ein- 
stein's eqs. has been presented here as a possible stable 
alternative to black holes for the endpoint of gravitational 
collapse. Realizing this alternative requires that a quan- 
tum gravitational vacuum phase transition intervene be- 
fore the classical event horizon can form. Although only 
the static spherically symmetric case has been consid- 
ered, it is clear on physical grounds that axisymmetric 
rotating solutions should exist as well. Since the en- 
tropy of these objects is some 20 orders of magnitude 
smaller than that of a typical stellar progenitor, there is 
no entropy paradox and instead a violent process of en- 
tropy shedding, as in a supernova, is needed to produce a 
cold GBEC or 'grava(c)star' remnant. The shell with its 
maximally stiff eq. of state p = pc 2 , where the speed of 
light is equal to the speed of sound could be expected to 
produce outgoing shock fronts when struck. These may 
serve to distinguish gravastars from black holes obser- 
vationally, and possibly provide a more efficient central 
engine for energetic astrophysical sources. The spectra 
of gravitational radiation from a struck gravastar should 
bear the imprint of its fundamental frequencies of vi- 
bration. Finally, we note that the interior de Sitter re- 
gion with p = —pc 2 may be interpreted also as a cos- 
mological spacetime, with the horizon of the expand- 
ing universe replaced by a quantum phase interface. 
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